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Spectroscopic probes of isolated nonequilibrium quantum matter: 
Quantum quenches, Floquet states, and distribution functions 
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We investigate radio-frequency (rf) spectroscopy, metal-to-superconductor tunneling, and angle- 
resolved photoemission spectroscopy (ARPES) as probes of isolated out-of-equilibrium quantum 
systems, and examine the crucial role played by the nonequilibrium distribution function. As an 
example, we focus on the induced topological time-periodic (Floquet) phase in a two-dimensional 
p + ip superfluid, following an instantaneous quench of the coupling strength. The post-quench 
Cooper pairs occupy a linear combination of “ground” and “excited” Floquet states, with coefficients 
determined by the distribution function. While the Floquet bandstructure exhibits a single avoided 
crossing relative to the equilibrium case, the distribution function shows a population inversion of the 
Floquet bands at low energies. For a realization in ultracold atoms, these two features compensate, 
producing a bulk average rf signal that is well-captured by a quasi-equilibrium approximation. In 
particular, the rf spectrum shows a robust gap. The single crossing occurs because the quench- 
induced Floquet phase belongs to a particular class of soliton dynamics for the BCS equation. 
The population inversion is a consequence of this, and ensures the conservation of the pseudospin 
winding number. As a comparison, we compute the rf signal when only the lower Floquet band 
is occupied; in this case, the gap disappears for strong quenches. The tunneling signal in a solid 
state realization is ignorant of the distribution function, and can show wildly different behaviors. 
We also examine rf, tunneling, and ARPES for weak quenches, such that the resulting topological 
steady-state is characterized by a constant nonequilibrium order parameter. In a system with a 
boundary, tunneling reveals the Majorana edge states. However, the local rf signal due to the edge 
states is suppressed by a factor of the inverse system size, and is spatially deconfined throughout 
the bulk of the sample. 
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I. INTRODUCTION namics can be induced in low temperature superconduc¬ 

tors on ultrashort timescales. Our results apply in a tran- 
A. Idea and overview of results sient window over which the time-periodic phase can be 

stabilized, before pair-breaking or other processes destroy 


In recent years, there have been numerous studies on 
Floquet systems as they possess interesting topological 
phases absent in static systems 00 . Most of these are 
induced by external driving. For example, a Floquet 
topological insulator may arise when graphene is irra¬ 
diated by circularly polarized light 0 i 1 S, 0- How¬ 
ever, a topological Floquet phase can also be created in 
an isolated system such as a BCS superfluid following 
an instantaneous quench of the interaction strength [8|- 
0. The periodic modulation of the order parameter is 
self-generated by dynamics, and arises due to the single 
unstable mode associated to the BCS instability of the 
normal state |UI IT~if . 

In either of these two cases, the Floquet (quasienergy) 
band structure and the distribution function (nonequilib¬ 
rium occupation number) are equally important. Much 
of the previous work focused on the connection between 
the band structure and topological properties 00 , but 
did not take into account the distribution function which 
might be essential to determine experimental observ¬ 
ables. The problem of irradiated graphene has been 
considered mainly in the Landauer-Buttiker formalism, 
wherein the mode occupation is fixed by ideal leads 
1 i I 0- Even in that simplified setting, however, 
the relationship between the topology and Floquet edge 
states to measured transport coefficients is complicated, 
and the latter are not generically quantized p,[S,0. More 
recent works have considered Floquet systems that are 
isolated or weakly coupled to the environment [IM3. 
If external driving is applied to a system initially in its 
ground state, the particles will populate several Floquet 
bands, leading to a nontrivial distribution function. 

To better understand the connection between the dis¬ 
tribution function and observables, in this work we study 
the radio-frequency (rf) spectrum [18j] of an isolated topo¬ 
logical Floquet system. We consider the topological p+ip 
Floquet superfluid induced by an instantaneous quench of 
interaction strength in a 2D system of spinless fermions 
0, @]. The long-time asymptotic steady state occupies 
two Floquet bands related by particle-hole symmetry. 
We determine the distribution function as well as the ex¬ 
plicit BCS dynamics of the steady state, and derive the 
expression for the rf spectrum. This out-of-equilibrium 
system could be realized in an ultracold gas of fermionic 
atoms interacting through a p -wave Feshbach resonance 


Assuming it can also be realized in solid, we discuss 
the possibility of probing this nonequilibrium system by 
normal metal to superconductor tunneling [19(, whose 
setup closely parallels that of rf. Finally we consider 
momentum-resolved spectra, as could be measured us¬ 
ing angle-resolved photoemission spectroscopy (ARPES) 
[20j . Recent THz pump-probe experiments [2lL [22l | have 
shown that order parameter (“Higgs mode”) quench dy- 


it [alii. 

We show that the bulk rf and tunneling spectra are in 
general very different, and this is tied to the nonequilib¬ 
rium distribution function. The rf and tunneling spectra 
are both one-particle observables, but there is a crucial 
distinction |23[. In an ultracold gas, rf radiation can in¬ 
duce internal transitions in the atoms between states that 
participate in pairing and states that do not. The rf sig¬ 
nal depends upon the occupation of both the paired and 
unpaired states. By contrast, tunneling from an idealized 
metal tip gives a current that is sensitive to the distribu¬ 
tion function in the superconductor only on energy scales 
of order £p , the Fermi energy in the tip. The tunnel¬ 
ing conductance at small bias is completely determined 
by the retarded single particle Green’s function in the 
Floquet phase. The same quantity determines transport 
in Floquet scattering theory 0, and signals the absence 
or presence of edge states in a system with a boundary 
0, M HH; however, it does not encode the occupation of 
bulk or edge states. 

In our quenched-induced topological superfluid, we 
show that there is exactly one avoided crossing in the 
Floquet bandstructure. At the same time, there is a 
population inversion in the occupation of these bands, 
so that the “upper” (“lower”) quasienergy band is occu¬ 
pied at low (high) energies. These two features compen¬ 
sate in the bulk rf prediction, so that the average sig¬ 
nal is well-approximated by that of a quasi-equilibrium 
p + ip superfluid. In particular, the spectrum exhibits 
a robust gap. Deviations from a quasi-equilibrium pic¬ 
ture appear as satellite coherence peaks, which have been 
studied previously jlS], To compare, we also consider the 
rf signal that obtains by populating only the lower Flo¬ 
quet band. In this case the gap disappears for sufficiently 
strong quenches that induce large variations in the order 
parameter (i.e., strong driving). The tunneling signal 
is ignorant of the distribution function, and also depends 
sensitively on the quench. For strong quenches there is no 
gap in the bulk tunneling spectrum. A quench-induced 
population inversion in a single band Floquet system was 
studied in [2(1. 


Unlike toy models sometimes considered in the litera¬ 
ture, the drive frequency in our quench-induced Floquet 
system is much smaller than the bandwidth. This im¬ 
plies that there is significant folding of the unperturbed 
spectrum across the quasienergy zone. Although this sit¬ 
uation should be the norm rather than the exception ex¬ 
perimentally i, it is usually more complicated to un¬ 
derstand than the opposite, high frequency limit 00. 
In particular, one typically expects small gaps to open 
every time the spectrum crosses the zone, and this can 
complicate predictions for experiment [1]. 

By contrast, the presence of a single avoided band 
crossing in our Floquet bandstructure is a consequence 
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of the BCS dynamics that generate it. In particular, the 
crossing arises because the quench-induced Floquet phase 
belongs to a particular class of soliton dynamics for the 
integrable BCS equations 0111. The population in¬ 
version is then a topological consequence of this, since 
the texture of the instantaneous Anderson pseudospin 
description of the BCS state is conserved by the dynam¬ 
ics d,Hj. In other isolated topological Floquet systems, a 
population inversion is also expected for an odd number 
of avoided crossings. 

We also consider the rf, tunneling, and ARPES signals 
that result from weaker quenches in the p + ip superfluid, 
such that the order parameter asymptotes to a nonequi¬ 
librium constant value. In addition to bulk spectra, we 
consider a semi-infinite geometry with an edge. We show 
that while tunneling reveals the Majorana edge states, 
local rf spectroscopy is ill-suited to find them. This is 
because the signal is suppressed by a factor of the sys¬ 
tem size, and is not spatially localized to the edge of the 
cloud. 


B. Review: Quench-induced topological Floquet 
phase in a p+ip superfluid 

In this subsection, we briefly review the quench in¬ 
duced nonequilibrium state 0 [§] whose rf spectrum is 
investigated. The system is governed by the BCS Hamil¬ 
tonian, 

7 2 

h =J2 W c k c k - 20 J2 k ' < l4 c -k c -q c q> ( L1 ) 

k k.q 

where Ck annihilates a fermion with momentum k = 
{k x ,k v } and mass m = 1. The prime in Y2' indicates 
that the summation runs over momenta with nonnega¬ 
tive y-components, i.e. k v > 0. G > 0 is the interaction 
strength. 

This Hamiltonian [Eq. (11.11) ] has a p + ip ground state 
whose order parameter assumes the form 

/ 

A(k) = -2 G ^2 k ■ q( c -q c q) 

q 

= A 0 (k x -ik y ) = A 0 ke-^ k , ( 1 . 2 ) 

where </>*, is the polar angle of 2D vector k. Ao is the 
order parameter amplitude. The quasiparticle energy of 
such a p + ip paired state is 

E k = -pf + k^A^ ( 1 . 3 ) 

Its spectrum is fully gapped as long as the chemical po¬ 
tential p ^ 0. The critical point p = 0 separates two dis¬ 
tinct topological phases —the weak pairing BCS phase 
(p > 0) and the strong pairing BEC phase (p < 0). The 
BCS phase is topologically nontrivial, i.e. if the system 
in this phase possesses a boundary, the Majorana edge 


state would appear. The squared order parameter (Ao) 2 
and chemical potential p carry units of density; the fixed 
particle density n sets the natural scale. All dimension¬ 
ful quantities for a quench can be expressed in units of n, 
with at most logarithmic dependence upon an ultraviolet 
energy cutoff A [9]. 

The system is initially prepared in the p + ip ground 
state of the pre-quench Hamiltonian with interaction 
strength Gi. Then the interaction strength is suddenly 
quenched to a different value Gf. After the quench, the 
system evolves as a superposition of many body eigen¬ 
states of the post-quench Hamiltonian and acquires a 
steady state as t — > oo. 

In general, the asymptotic long-time evolution of the 
order parameter amplitude can be expressed as 

A{t) = A 00 {t)e~ 2i ^\ (1.4) 

where Aoo(t) could be 0 (phase I), a nonzero constant 
(phase II) or an oscillating function with time period T 
(phase III), depending on the strength and direction of 
the quench 0 9i, 0, 0 • In phase II, the real constant 
p ^ plays the role of a nonequilibrium chemical potential 
and determines the topological properties of the quench 
induced state. It is topologically nontrivial (trivial) when 
Moo > 0 ( p oo < 0). The out-of-equilibrium phase diagram 
is shown in Fig. [T] where we use A^' 1 and the order 
parameter amplitudes of pre- and post-quench Hamilto¬ 
nians’ ground states, to indicate the coordinates of the 
quench @, |i|. 

The main interest in this work is phase III, which can 
be understood as a quench-induced Floquet topological 
superfluid phase Q. For both phases II and III, we use 
the term “topological” to denote the presence or absence 
of Majorana edge modes in a sample with a boundary, 
as encoded in the single particle retarded Green’s func¬ 
tion 010. This notion of topology relates to the 
spectrum of excitations that can be reached from the 
quench-induced nonequilibrium state by the application 
of a weak spectral probe, i.e. a frequency-dependent mea¬ 
surement. There is a different notion of topology that de¬ 
scribes the many-body wavefunction itself, which is here 
always a BCS state. This is the pseudospin winding num¬ 
ber that can be observed by an equal-time measurement, 
such as time-of-flight, and which cannot change follow¬ 
ing a quench. The conserved pseudospin winding does 
not determine whether or not edge states will appear; it 
instead encodes the occupation of the states ii- Invari¬ 
ance of the winding number of the state was also noted in 
externally driven Floquet 0 and quenched topological 
0 models. 

In the following, we focus on weak-to-strong quenches 
within phase III where the order parameter amplitude 
is given by a Jacobi elliptic function The constant 
Aqq in phase II can be considered as a trivial case with 
vanishingly small oscillation amplitude. 
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by the BCS product state 


FIG. 1: (color online). Quench phase diagram |3l[ show¬ 
ing 3 distinct dynamical phases characterizing the out-of¬ 
equilibrium states induced by sudden quench of the interac¬ 
tion strength in a 2D p + ip fermion superfluid US]. The 
order parameter amplitude A oc (t) either vanishes, converges 
to a nonzero constant, or oscillates persistently for quenches 
in phase I, II, or III, respectively. The vertical and horizon¬ 
tal axes show Aq 1) and A ( /\ the order parameter amplitudes 
associated to the ground states of the pre- and post-quench 
Hamiltonians. Aqcp indicates the critical point where the 
ground state chemical potential vanishes. It separates the 
topologically nontrivial BCS and trivial BEC phases in equi¬ 
librium. The purple dashed curve is its nonequilibrium exten¬ 
sion. Each point to the left (right) of this line in II represents 
a topologically nontrivial (trivial) nonequilibrium state with 
positive (negative) . Phase III is a quench-induced Floquet 
topological superfluid state, which is our main focus here. The 
labeled points are particular quenches discussed in the text. 
The amplitudes A and Aqcp are measured in units of the 
inverse length y/n, where n is the fixed particle density [gj. 


C. Floquet states and distribution function 


f 

I^W) = n [«k(i) + Wk(i)4 c - k 1°). ( L6 ) 

k 

where the time-dependent coherence factors follow the 
Bogoliubov-de Gennes (BdG) equation 


dt 


Uk(t) 
Vk (t) 


2 

kA a 


-/T* 

(t) 


kA* 

kL 1 

2 


,(*) 

Moo 


u k (t) 
Vk (t) 


(1.7) 


When A 00 (t) is periodic in time as in phase III @ with 
A 00 (t + T) = A^jt), Eq. (11.711 has a pair of Floquet so¬ 
lutions. These are states taking the form |\|/( F )(t)) = 
I & F \t))e lE( )t , where shares the same period 

T with A oc (t), + T)) = |3>^(i)). The quasi¬ 

energy E( p ) is defined up to integer multiples of the fre¬ 
quency Q = 2‘k/T. A Floquet state is a solution subject 
to the special initial condition that |^( F )(0)) is an eigen¬ 
state of time evolution operator U(T ) over one period: 
U(T)\^ f \ 0)} = e iE(F)T \^^ F \0)). 

We will always assume a particular type of physical 
initial condition for the many-fermion system, which is 
the p + ip ground state of the pre-quench Hamiltonian. 
The coherence factors in the asymptotic steady state are 
then given by a superposition of two Floquet states, both 
of which solve the BdG equation Eq. m, that is 


Ukit) 

/i-Tk 

g 

C-+. 

_1 

Vk (t) 

V 2 

1 

1_ 


a +iEi F h 


./ 1 + 7k 


V 2 

r ul {F \t)_ 


( 1 . 8 ) 


c ,-iE^ ) t+iT k 


These two Floquet states are related by particle-hole 
symmetry, and the relative coefficients are set by the 
“Cooper pair distribution function” 7k (—1 < 7k < 1), 
which determines the nonequilibrium occupation number 
/ k QP) of the Floquet bands p, HH . 


These are related via 


The order parameter amplitude Ajf) is determined 
bv the self-consistent BCS dynamics m 0 ns m 
1301 ]. The dynamics of individual Cooper pairs or (pairs of 
quasiparticle excitations) can then be obtained by solving 
the effective mean field Hamiltonian, 

//B J C = El^“^ (4Ck + ClkC " k) i 

k [ + fcA„(i)Ckcl k + ftA^(t)c_ k CkJ 

(1.5) 


fr = l( l+7k), (1-9) 

so that 7 k = —1 (+1) corresponds to the absence (pres¬ 
ence) of a pair of excited quasiparticles occupying states 
±k. This is different from the occupation of the bare 
fermions ((c k c k + c^ k c_ k )). 

We will demonstrate in the following sections that the 
distribution function possesses topological information of 
the pre- and post- quench states. It can be extracted from 
the rf spectroscopy but not from the superconductor- 
normal metal tunneling experiment. 


Here we have absorbed the polar phase of A(k) into the 
fermion operators via c_ k c k —> C- k c k- This 

transformation also boosts into the “rotating frame,” so 
that /To,, is transferred from the order parameter ampli¬ 
tude to the kinetic term. 

The corresponding many body wavefunction is given 


D. Outline 

This paper is organized as follows. In Sec. im we de¬ 
rive general expressions for rf spectroscopy, tunneling am¬ 
plitudes, and ARPES employed to probe the quench- 
induced out-of-equilibrium system. In Sec. cm as a 
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warm-up we compute the rf, tunneling, and ARPES spec¬ 
tra associated with several quenches in phase II with var¬ 
ious topological properties. We examine the possibilities 
of extracting the topological character of both pre- and 
post- quench states by analyzing the rf spectrum. In 
Sec. nn we turn to the analysis of the Floquet phase 
III. We study thoroughly two representative phase III 
quenches: one located in the vicinity of the phase II III 
boundary, while the other lies deep in the phase III with 
weak initial interaction strength. The former (latter) is 
characterized by small harmonic (large anharmonic) am¬ 
plitude modulations in A oc (t). We present and discuss 
the results for the rf, tunneling, and ARPES spectra, 
and we also discuss the connections between the Floquet 
bandstructure and the integrable BCS dynamics. Finally, 
in Sec. [V] we compute the local rf signal of a phase II post- 
quench system with a boundary, looking for signatures of 
Majorana edge states. We give our conclusion in Sec.lVIl 
Expressions for rf and tunneling harmonics relevant 
to probing the Floquet phase are relegated to Appendix 
El The detailed derivations of the explicit Floquet state 
wavefunctions as well as the distribution function associ¬ 
ated with phase III quenches are given in Appendix [B] 


II. RF, TUNNELING SPECTROSCOPY, AND 
ARPES 


A. Hamiltonian and current 


In an rf experiment, we assume that the system un¬ 
der study is realized in an ultracold fermion gas that 
possesses two relevant hyperfine states: |1) which partic¬ 
ipates in pairing, and |2) which does not. 

The Hamiltonian for atoms in the non-pairing state |2) 
is given by 


H, 


(d) 


y. (v+£2,i)4.d k , 


( 2 . 1 ) 


of 12) atoms at position rp, 


n d { r 0 ) = d t (r 0 )d(r 0 ) = ^ [d kl dk 2 e* (kl k2)r ° 

ki,k 2 


(2.3) 


while N d is the total number of atoms in state |2) 

Ad = ^<i k dk- (2.4) 

k 


Aside from the rf, we also consider normal metal to 
superconductor tunneling, assuming that the same post- 
quench asymptotic steady state can be hypothetically re¬ 
alized in a superconductor. The quenched condensate 
is brought into contact with a normal metal tip, whose 
Hamiltonian is given by 

H i o d) = J2{ ! T + V ) d i d ^ (2.5) 

k 


Here the tunneling bias V is the energy of the tip rela¬ 
tive to the condensate. We assume that the tip is ideal, 
characterized by a constant density of states uq, and that 
prior to contact only states with momentum k < qp are 
occupied. The Hamiltonian 


H t =T [c t (r 0 )d(r 0 ) + d t (r 0 )c(r 0 )] 

= r J2 [4<v i(k_q)ro +<c k e- i(k - q)ro 

k,q L 


( 2 . 6 ) 


describes the tunneling between the metal tip and the 
superconductor at the contact point r 0 . The observable 
is again the local current at this point [Eq. (12.31) ]. 

Comparing the Hamiltonian of tunneling [Eq. (12.61) ] 
with that of the rf coupling [Eq. (12.2|l ] in the rotating 
frame d k —> d k e _ *“ z ' t , the only difference is that while 
the contact is global for rf, it is local for tunneling. They 
are identical if the transfer matrix element T in Eq. (12.61) 
is replaced with TA k . q , meaning the momentum is con¬ 
served in an rf experiment but not in tunneling. 


where £ 2,1 denotes the atomic transition energy between 
the states |2) and |1). This energy separation is typically 
much larger than any relevant to the many-body dynam¬ 
ics. d k annihilates an atom in state |2) with momentum 
k. In a p- wave paired superfluid, the rf radiation induces 
a transition in an atom with momentum k from state |1) 
to 12); an unpaired state |1) atom with momentum —k is 
left behind. This process is described by the transition 
(coupling) Hamiltonian 

Ht = TJ2 [e iULt cld* + e^dU] , (2 . 2) 

k 

where lot is the frequency of the rf field. 

We will consider two observables: the local rf cur¬ 
rent /(ro) = (drid (rg) / dt) and the global current I = 
(dNd/dt). Here nd(ro) is the spatial density distribution 


B. General rf and tunneling formulae 


Using linear response theory, we obtain the rf current 
as a function of the coherence factors 


I rf (t) =27^Re^4, q f dt'e^-^'-V 

k.q d -°° 




q u k\ 


(2.7) 


where to = lol — £ 2,1 is the detuning frequency, and is 
the occupation number of atoms in state |2) before per¬ 
turbed by rf radiation. The bare energy in the rotating 
frame is given by = k 2 /2 — /i„. 

Similarly we can obtain the expression for the tunnel¬ 
ing current by replacing to and d k q in Eq. m with — V 
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and 1, respectively. In addition, now = 9{qp — q ), 
where 9 denotes the Heaviside unit step function. The 
differential conductance defined as G{t) = — is 

given by 


G(t) = 2T 2 v 0 Re J2 [ dt’e 

x K( i ')«k(0 + Wk(t)Mk(i')] (2-8) 

= - 2TV 0 Im ^ f dt gw ( t) ^ 

k 


where V = — V — t; qF . On the last line, we’ve used the 
fact that 

K(*>kO) + ul(t)u k (t')} = ({ck(t),4(t')}) = 

is the single particle retarded Green’s function for a BCS 
product state. 

The rf current [Eq. (12.71) ] depends explicitly upon the 
distribution of the |2)-state atoms and, as we will see, 
that of the paired ones as well. The tunneling conduc¬ 
tance [Eq. (12.81) ] by contrast depends only on the general¬ 
ized spectral function (Wigner transform of the retarded 
Green’s function) in the condensate. 


C. RF and tunneling for a Floquet system 


To find the rf signal from a topological Floquet system 
induced by phase III quench, we insert Eq. m into 
Eq. (12.7II . and obtain the time-averaged rf current 


lr f = nT ‘ 2 Y. 

n. k 


< 

+ 


a-7k)(i-/^ ) )-(i+7k)/^]' 
x |un,k| 2 <i(w - £ k - El F) + nfl) 

"(l + 7k)(l-/^ ) )-(l-7k)/^ ) l 


, X |u n , k | 2 d(w - £ k + - nGl) , 

(2.9) 

where u n ,k and v Ui k are the Fourier coefficients of u ^ ; ( t ) 
and 


.(F) 

,(F) 


(*) 

W 


E 


^n,k 

^n,k 


g —inflt 


( 2 . 10 ) 


In a Floquet phase, the rf current and tunneling con¬ 
ductance also show periodic modulation at harmonics of 
the drive frequency f1. In Appendix [A] we compute 

Irfip) = £ §e* m Irf(t )> G(p) = £ ^e ip ™G(t). 

( 2 . 12 ) 

These harmonic amplitudes provide additional observ¬ 
ables that can be used to characterize a Floquet phase. 

The time-averaged rf amplitude in Eq. (12.91) depends 
explicitly on the distribution function 7 k, irrespective of 
the |2)-state occupation The same conclusion ap¬ 

plies for all harmonics [Eq. (I A 11) ]. By comparison, even 
the instantaneous tunneling conductance G(t) is indepen¬ 
dent of 7k, if the tunneling probe has a constant density 
of states. 


D. ARPES formulae 


In ARPES experiments, one measures the Wigner 
transform of the single particle lesser Green’s function 
Gk,<(t,t') 


S{u, Mo) = J dte + , 

(2.13) 

which equals 

-iGk,<(t,t') = (cl(t')c k (t)) = v£(t')v k (t). (2.14) 


Inserting Eq. (12.101) into Eq. (12.131) and Eq. (12.141) . we 
arrive at the time-averaged ARPES signal from a quench 
induced Floquet system, 


SV,k) 



(1 - 7k)|un,k| 2 ^(-w - eI F) + nD)! 

+ (1 + 7k) |!ln,k| 2 (f( iO + E £ ^ — ntt)j 

(2.15) 


III. WEAK QUENCHES WITH CONSTANT 
ORDER PARAMETER: BULK SPECTROSCOPY 
IN PHASE II 


A. Phase II introduction 


A version of Eq. (12.91) assuming f^ = 0 appeared previ¬ 
ously in [ij| . Here all the oscillating terms are discarded 
since their time averages vanish. From an analogous cal¬ 
culation, we determine the time-averaged differential con¬ 
ductance, 

G = 27 t 7 ^„EE[ |5 ”’ k|2 ^“^ +nn) . 

nk L+Kk| 2 6 (V + E ( k F) - 7 ^n) 

( 2 . 11 ) 


For the long time asymptotic steady state following 
a phase II quench, only the zeroth order Fourier coeffi¬ 
cients uo.k,'5o,k are nonzero. The time-averaged rf cur¬ 
rent, tunneling conductance, and ARPES signal can be 
easily determined from Eqs. (12.91) . (12.111) . and (12.151) . us¬ 
ing the exact results for A„, , and y(fc) that can be 

obtained for a particular quench |9|. In what follows, un¬ 
less otherwise stated, the non-pairing state is assumed to 
be initially unoccupied for all momenta in the rf calcula¬ 
tion, fl d) = 0. 
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FIG. 2: (color online). Bulk rf spectra (red solid curves) 
for nonequilibrium states induced by (a) BCS to BCS quench 
“a,” (b) BEC to BCS quench “b,” as indicated in the quench 
phase diagram Fig. [l] The green dashed curves show the rel¬ 
ative weight imparted by the distribution function equal to 
1 , 1+ry ( kuJ ' > in the positive and negative frequency do¬ 

main, respectively. The distribution function yt, is illustrated 
in the inset. It winds from +1 at k = 0 to -1 as k — > 00 
for the BEC to BCS quench b, but goes from -1 to -1 for 
BCS to BCS quench a. For both quenches, the gap is located 
between — and 0 in the spectrum, indicating the (topolog¬ 
ically nontrivial) weak-pairing BCS nature of the post-quench 
state. A discontinuous jump appears at u> = 0 in (a) which is 
associated to a quench from a BCS initial state. The coordi¬ 
nates for quenches a and b are {A^, Ajy^} = {0.8, 0.5}Aqcp 
and {1.1, 0.8}Aqcp, respectively. 


FIG. 3: (color online). Bulk rf spectra for nonequilibrium 
states induced by (a) BEC to BEC quench “c,” (b) BCS to 
BEC quench “d,” as indicated in the quench phase diagram 
Fig CD As in Fig. [2] the red solid and green dashed curves 
respectively correspond to the actual spectrum and relative 
weight, while the inset shows the distribution function 7 k- 
The distribution function contains an even and odd number 
of zeros for quenches c and d, respectively. The post-quench 
states corresponding to both (a), (b) are topologically trivial, 
i.e. there would be no Majorana edge modes in a system with a 
boundary, and the non-equilibrium phase is “BEC-like.” This 
is indicated by the gaps which appear between 0 to 12/tool. In 
(b), the rf spectrum exhibits a discontinuous jump at w = 
0, indicative of the BCS initial condition. The coordinates 
for quenches c and d are {A^ 1 , Aj/ 5 } = {1.02, 1.2}Aqcp and 
{0.7,1.1} Aqcp, respectively. 


B. Phase II: RF 


The rf spectrum of post-quench state with constant 
order parameter takes a nonzero value only within the 
region 


w e 


[—2/too, — U [0, 00 ), 


A 



w e [-Ai,-2/rJ U [0,oo), -^>AGo>0, 
we [—Ai,0]U[|2 Moo |,oo), /too < 0. 


(3.1) 


It is composed of two continuous parts separated by a 
gap of width 12/too | or A^. The gap appears for negative 
(positive) w for topologically non-trivial BCS-like (trivial 
BEC-like) states with /too > 0 (/<co < 0) [1, [24[ . This 
feature in principle allows one to distinguish topological 
and trivial phases via a bulk measurement, even in the 
ground state pi} . 

A crucial difference from tunneling is that the weight 


always extends to zero on one side of the rf gap. This is 
due to the fact that fermions with k —> 0 are not coupled 
to the order parameter A for p- wave pairing, and thus 
transition between hyperfine states at the bare frequency. 

The component on the positive (negative) detuning fre¬ 
quency side is due to the process where an rf photon 
with energy u>l breaks the ground (excited) Cooper pair 
labeled by k^ and excites one of the atoms to the dif¬ 
ferent internal state 12 ), which does not participate in 
pairing. An unpaired state |1) atom is left behind j§}. 

Here k u = the solution to w = ±£4 + 

which comes from the conservation of energy. The term 
describing this process is weighted by 1 

which reflects the occupation number of the ground (ex¬ 
cited) Cooper pair. 

As required by the conservation of pseudospin winding 
number Q = s 2 (0) — s z (oo), the distribution function 
y(fc) must wind from +1 at k = 0 to —1 at k —> 00 if the 
pre- and post-quench states are in different topological 


1+7(M 
2 ’ 
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FIG. 4: The bulk tunneling spectra of the nonequilibrium 
state following quenches (a) “a,” (b) “f,” (c) “d,” as indicated 
in the quench phase diagram Fig. [I] The corresponding coor¬ 
dinates are {Aq 1) , A]/'} = {0.8, 0.5}Aqcp, {0.8, 0.9}Aqcp, and 
{0.7, 1.1}Aqcp, respectively. In (a), the long time asymptotic 
nonequilibrium state has > A^; two coherence peaks ap¬ 
pear at V = iVmin - A^, > /Too > 0 and /i„ < 0 for quenches 
f and d respectively. In (b), (c), the tunneling signal jumps 
from 0 at one edge of the gap V = and grows continuously 
at the other edge V = —/Too- The discontinuous jump occurs 
at the right (left) side of the gap in Fig(b) [Fig(c)]. 


FIG. 5: (color online). The ARPES spectra of the nonequi¬ 
librium state induced by quenches (a) “a,” (b) “b,” (c) “c,” 
(d) “d.” 



0.03 


0.02 


0.01 


0 


phases, i.e. if Mo< 0 Q. When MoVoo > 0, 7(k) 
does not wind but approaches — 1 as k —> {0, 00}. In 
particular, the difference between these two cases is the 
distribution function value at k = 0 : 


limiO) 

k —>0 


- 1 , > 0 , 

+1, ztfV- < 0. 


(3.2) 


The rf signal from which the distribution function can 
be extracted therefore possesses topological information 
of the state before and after the quench. A discontinuous 
jump at uo = 0 would appear in the rf spectrum if and 
only if /zj 75 > 0. This discontinuous jump exists at the 
right (left) side of gap edge if Moo > 0 (m=o < 0). On the 
other side of the gap, the spectrum grows continuously 
from 0. This introduces the feasibility of using quantum 
quench in the detection of topological information of the 
system. 

As an example, we evaluated the time-averaged rf spec¬ 
tra of the steady states induced by four different phase 
II quenches. These quenches are indicated as points “a” 
(BCS to BCS), “b” (BEC to BCS), “c” (BEC to BEC) 
and “d” (BCS to BEC) in Fig. [T] For quenches a and b 
(c and d), the asymptotic chemical potential > A^/2 
(/Zoo < 0), and the corresponding rf spectra are illustrated 
in Fig. [2] (Fig. [3]). The gap spans from — to 0 in Fig.[2l 
and from 0 to —2/4^ in Fig. [3l as expected. In addition, 
the discontinuous jump appears at co = 0 for quenches a 
and d, which obtain from a BCS initial state. 


C. Phase II: Tunneling 

As in rf, a gap appears in the tunneling spectrum. 
Different from rf, it spans the symmetric interval from 
— E m \ n < V < E m - m , where E min is the minimum excita¬ 


tion energy (spectrum gap): 


[ ^min 1 

Moo 

> Ai, 


-^min — \ /^oo 1 

Al 

> Moo > 0, 

(3.3) 

{ -M=, 

Moo 

<0, 


Vmin = A^-y/2 /Z c 

0 - 

AF 

(3.4) 


Furthermore, the topological properties of the post- 
quench state can also be inferred from tunneling sig¬ 
nal near the gap edge V = ±E m \ n . When /z„ > A^, 
G(V) exhibits two coherence peaks around V = iV^n, 
as shown in Fig. 0] (a). The coherence peaks disappear 
for /z x < A^, and instead the spectrum shows a dis¬ 
continuity on the positive (negative) edge of the gap for 
Hao > 0 (BCS-like) [/z^ < 0 (BEC-like)], see Figs. |4{b) 
an dStc). However, unlike the rf, the tunneling spectrum 
does not reveal information about the pre-quench state 
since the distribution function disappears from the dif¬ 
ferential conductance. 


D. Phase II: ARPES 


There are two branches in the ARPES spectrum of 
the post-quench state with constant order parameter. 
The upper (lower) one results from the excited (ground) 


1—r(fc^) 

2 


For 


Cooper pair, and is weighted by 1+7 ^ fc ^ 

a ground state, the upper branch disappears. Simi¬ 
lar to tunneling, there is no signal within the region 
co £ {—E m i n , E m i n ) [see Eq. (13.31) ]. From Eq. (12.151) . the 
intensity of the signal is determined by both the distribu¬ 
tion function 7k and the moduli of the coherence factors 
|?zo,k|) |uo,k|- Therefore, the distribution function is also 
measurable from an ARPES experiment. 

In Fig. [5l we plot the ARPES signal for the steady state 
following phase II quenches a, b, c, and d. To compute 
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the signal, we replace the delta function in Eq. (12.151) 
by a smearing function S sm (x) = n ^ x 2 +rj 2 ^ , and set 77 = 
(6 x 10“ 4 )£f, where £f = 27rn is the Fermi energy of 
the system (n is the density). For quenches a and b, the 
post quench state is topologically nontrivial > 0 ) 
and, as a result, |f(o,k| = 0 at k = 0. Examining the 
ARPES spectrum at k = 0, we find it reaches maximum 
at u) = —/Koo for BCS to BCS quench a [Fig [SJa)] and 
disappears for BEC to BCS quench b [Fig [5j b)]. The 
absence of signal for quench b at k = 0 is due to the 
population inversion that occurs when the pre- and post- 
quench state are in different topological phases. On the 
other hand, the post quench states induced by quenches 
c and d are topologically trivial < 0 ), and |uo,k| = 0 
at k = 0. The ARPES signal with k = 0 maximizes 
at w = —/Too for BCS to BEC quench d [Fig[5][d)] and 
vanishes for BEC to BEC quench c [Fig[5][c)]. For quench 
d, the peak at to = — is a signature of the population 
inversion. 

Based on the discussion above, we conclude that topo¬ 
logical information of the state before and after the 
quench can be inferred from the ARPES signal. The 
pre- quench state is in BCS (BEC) phase when the sig¬ 
nal is present (absent) at k = 0. In addition, if the peak 
occurs at negative (positive) frequency, the post-quench 
state is topologically nontrivial (trivial). 



t/T 


IV. STRONG QUENCHES WITH PERIODIC 
ORDER PARAMETER: BULK SPECTROSCOPY 
IN FLOQUET PHASE III 


A. Phase III introduction 


The problem of computing the long-time evolution of 
a paired fermion superfluid following a quantum quench 
can be reduced to solving the integrable dynamics of 
coupled Anderson pseudospins. Following a similar con¬ 
struction employed in the study of s-wave superfluids 
HHim, 113,0, 0], in a Lax spectral method is used to 
find the analytical solution for the 2D p+ip model stud¬ 
ied here. For a system of N pseudospins, one introduces 
the spectral polynomial Q 2 n(u), a conserved integral of 
motion for any value of it. In the limit N —> 00 , Q 2 n{u) 
respectively exhibits zero, one or two isolated pairs of 
roots in phases I, II, III of the quench phase diagram, 
Fig. CO Isolated roots are those separated from the posi¬ 
tive real axis in the complex u plane. These always come 
in complex conjugate or negative real pairs. The remain¬ 
ing roots give rise to a branch cut along the positive real 
axis in the thermodynamic limit. 

The isolated roots in turn parameterize the dynam¬ 
ics in an effective “reduced” problem of zero, one, or 
two collective pseudospins. The asymptotic evolution of 
the order parameter reaches a steady state described by 
Eq. (II.41) . and this is completely determined by this re¬ 
duced solution. The dynamics of individual pseudospins 
can then be reconstructed following a “bootstrap” pro¬ 
cedure, solving the BdG Eq. (E3 and exploiting the 
conservation of Q 2 n(u) to determine the occupation of 
the non-equilibrium spectrum by ground- or excited-state 
pairs. The explicit solution for phases I and II was pro¬ 
vided in [§] . The analogous calculation for phase III gives 
the combination of Floquet states in Eq. dm); the deriva¬ 
tion is relegated to Appendix iBl 

For most quenches in phase III, the spectral polynomial 
possesses two complex conjugate pairs of isolated roots: 
one pair ui,± = u\ tX ±iu\^ with iip r > 0 exists only in this 
phase, while the other it 2 ,± = ii 2 jr ± 1112,1 remains isolated 
in both phases II and III. Here u a}X and u a ,i indicate the 
real and imaginary parts of u a .±, a = 1,2. We define 






a = 1 , 2 , 


(4.1) 


FIG. 6: (color online). Oscillating order parameter A 00 (t) 
induced by two phase III quenches, labeled “A” and “B” 
in the quench phase diagram Fig. |T] The left panels show 
the time dependence of the modulus and phase, while the 
right depicts the A 00 (t) orbits in the complex A plane. 
The magnitude of lA^t)! oscillates between the two turn¬ 
ing points A 2 ± Ai [c.f. Eq. m\. a 2 is indicated by the 
center point of each orbit on the right. The coordinates for 
quenches A and B are {Aq 1> , A(/ ; } = {0.119, 0.514}Aqcp, and 
{0.00651, 0.825}Aqcp, respectively. 


with \u a \ = yWo.r + ^ai the modulus of root it a ,±. 

A 00 (t) = A oc (t + T) can be expressed in terms of the 
Jacobi elliptic function cn(z\M) and the four parameters 
A 12 and fi\ t 2 - The latter are functions of the phase III 
quench coordinates {A™, A)/ 4 }, the fermion density n, 
and the ultraviolet cutoff energy A. The period T is also 
determined by the isolated roots Si- In terms of the 
quench {A™, AJ/'} for A™ -C and A)/ 5 < Aqcp, it 
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is given by [§} 

8 ,n A ^ y 2 / 4«_ (A (/) )2 

A A<V2y , -(A«) 2 J ’ 

(4.2) 

where /x^ and denote the ground state chemical 
potentials corresponding to the pre- and post-quench 
Hamiltonians, respectively. In phase II, where 1 x 2 ,± is 
the only pair of isolated roots, /x^ = /x 2 and = A 2 , 
a constant. 

In the complex A plane, the orbit of A oc (t) encircles 
A 2 between the two turning points A 2 + Ai and A 2 — 
Ai. Fig. [6] shows the A«,(i) orbits of two different phase 
III quenches “A” and “B,” as indicated in the quench 
phase diagram Fig. [I] The relatively weak quench A 
is close to phase II, and it exhibits a small orbit, with 
harmonic time-dependence of the amplitude and phase. 
The other quench B is located deep in phase III. Since the 
oscillation amplitude is of the same order as the average 
value, this constitutes an example of “strong driving” 
with respect to the induced Floquet bandstructure. The 
amplitude and phase of the order parameter are strongly 
anharmonic in time for quench B. In the remainder of 
this section, we use these two points as representative 
quenches in phase III. 


T ' 


a$'V 2 /4 /) — 


B. Floquet bandstructure from BCS dynamics 


The explicit solution to the BCS dynamics of the Flo¬ 
quet system induced by a quench in phase III takes the 
form given by Eq. (11.81) . It is a superposition of two Flo¬ 
quet states related by particle-hole symmetry. We deter¬ 
mine the exact expressions for these Floquet states and 
the quasienergy spectrum in Appendix [B] We also show 
that the modulus of the distribution function | 7 (fc)| takes 
exactly the same form as previously obtained for phases 
I and II [§]. 

It will prove very useful to construct “phase II” 
(static quasiequilibrium) analogs of the Floquet states 
and quasienergies, which will serve as a reference point 
for comparison. We introduce 

u 2 (e) 

M e ). 

which correspond to the ground and excited Cooper pair 
states for a phase II quench with A^ = A 2 and /x^ = p, 2 . 
Here we define 


iE 2 {e)t 

u 2 (e) 


-u 2 (e) 


0 —iE 2 {e)t 


(4.3) 


E 2 (e) = \J (e/2 — /x 2 ) 2 + eA 2 , 

/ x /l e/2 - /x 2 

« 2 (e)=J- + ' / 


V2(e) = -1 


2 2 E 2 {e) ’ 

(l e/2 — /x 2 

2 2 E 2 {e) ' 


(4.4) 



E 



£ 


FIG. 7: (color online). Quasienergy spectrum E <1 '" 1 (e) 

in the extended zone scheme for quenches (a) “A” and (b) 
“B”, as indicated in the quench phase diagram Fig. [I] Here 
e = fc 2 , where k is the momentum. For both quenches, the 
quasienergy E^ f \e) (red solid curve) is quite close to the 
static phase II approximation for the quasiparticle excita¬ 
tion energy E 2 {s) (green dotted curve) for large £, and to 
—A 2 (e) + C (blue dashed curve) for small e. Here O denotes 
the oscillation frequency of A„ c (t), determined by the quench. 
The crossover of E^ (e) between “ground state” and “excited 
state” branches occurs around e±, defined in the text. 


In this section, we denote e = k 2 for modes with momen¬ 
tum ±k. 

In Fig. [3 we plot the quasienergy E^ f \e) (red solid 
curve) as a function of e = k 2 in the extended zone 
scheme for quenches A and B. The green dotted curves in 
these figures show the dispersion of the excited states in 
the phase II approximation E 2 (e), while the blue dashed 
curves are the ground state energies shifted by the os¬ 
cillation frequency —E 2 (e) + Q [Q = 2tt/T , where T is 
the quench-induced oscillation period of A^ (t)]. We find 
that in both cases the quasienergy spectrum exhibits a 
single “avoided crossing” or “Floquet bandgap,” in which 
the behavior of E^{e) crosses over from the ground to 
excited state phase II approximations when the separa¬ 
tion between —E 2 (e) + Q and E 2 (e) is small [Fig. [3- 
In particular, the quasienergy is close to E 2 {e) when 
e O £_, and to —E 2 (e) + Q when e -C £+. Here £± 
(e+ < £_) denote a pair of intermediate single particle 
energies determined by the quench [see Eq. (1B20I) for ex¬ 
plicit formulae]. Technically, £+ (£-) marks the energy at 
which the coherence factor w,l F )(£,t) [?/ F )(£,t)] reaches 
zero within one period [xx^ F ^(£ + , 0) = 0, i/ F )(e_, -j) = 0; 
see Fig. [3 discussed below]. 
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0 E + S k\ 

s 

FIG. 8 : (color online). Oscillation of the first component 
of the phase III Floquet solution associated to 

quenches A (a) and B (b). The top green (bottom orange) 
solid curve shows the maximum (minimum) of |M^ F ^(e,t)| 
within one period. The yellow (shaded) area enclosed by 
these two curves marks the region swept out by the periodic 
modulation. We compare these to the static phase II approx¬ 
imation for the coherence factors |« 2 (e)| (pink dotted curve) 
and 1(e)| (cyan dashed curve). The magnitude of the Flo¬ 
quet component |u (F ^(e, t)| oscillates within a narrow region 
around the phase II approximation for most single particle en¬ 
ergies e, and switches branches near e±, where the oscillation 
amplitude is maximized. 


The coherence factors of the Floquet state 
{w/ F )(e, t), v^ F \e, t)} oscillate little around those of the 
phase II approximation “ground state” {^(e), t >2 (e)} or 
“excited state” {^(e), —^(e)} BdG spinors for most 
single particle levels, switching branches in the same 
region where the quasienergy changes behavior. Fig. [ 8 ] 
depicts the oscillation of (er, t) | for quenches A and 

B. The yellow (shaded) area marks the region swept out 
by |it( F )(e,t)| in one period. It shows that |?/ F )(£, f)| 
undergoes small oscillations near |f 2 (e)| (cyan dashed 
curve) or |rt 2 (e)| (pink dotted curve) when e -C £+ or 
e>e_, respectively. 

As shown in Eq. (12.91) above, the Floquet bandstruc- 
ture itself contains only part of the information rele¬ 
vant for rf spectroscopy. The other ingredient is the 
non-equilibrium distribution function 7 (e), as appears 
in Eq. CLED- This sets the weights of the “ground” 
and “excited” Floquet state solutions proportional to 
exp [±i£'( F ) (e)t ]. The explicit formula for 7 (e) is identi¬ 
cal to phases I and II; the latter was obtained in I. In 
Appendix[B] we confirm that the combination of Floquet 


-1 —„- L_- 

0 £ + £- k) 

S 

FIG. 9: Distribution function 7 (e) for quenches (a) A and (b) 
B. Due to the conservation of the pseudospin winding number 
[ 9 ], both wind from +1 at e = 0 to —1 at e —> 00 . 

and occupation factors in Eq. (HD gives the correct fixed 
particle density n and self-consistent expression for A(t) 
in Eq. (11.41) : see Eq. (IB29I) . 

We plot the distribution functions for quenches A and 
B in Fig. [9] Similar to the quasienergy spectrum and 
coherence factors, 7(e) exhibits a crossover from —1 at 
large e to +1 as e —> 0. This holds true for any phase 
III quench, and implies that the non-equilibrium distribu¬ 
tion function shows a population inversion of the Floquet 
states at low energies. I.e., the “lower” Floquet band is 
occupied at large energies, but the “upper” one is filled 
for k < &f, c.f. Eq. (11.91) . 

In fact, the “winding” of 7(e) from minus one to plus 
one as e decreases from infinity is required by the topol¬ 
ogy. In particular, the pseudospin winding number that 
characterizes the instantaneous BCS state of the manv- 
fermion system cannot change following a quench [9|. 
(This is different from the question of Majorana edge 
modes, which can appear or disappear; these are en¬ 
coded in the retarded Green’s function winding num¬ 
ber sum, and this quantity can change in a quench 
across the quantum critical point Q|.) Given the single 
crossover of the quasienergy and Floquet coherence fac¬ 
tors relative to the static quasiequilibrium case [phase 
II approximation, Eqs. (Id. 3D and (HD], conservation of 
the pseudospin winding number for BCS initial states 
implies that 7(e) must go to +1 as e 0 , so that 
s k = ( c k c k + c _k c -k ~ l )/ 2 = + 1/2 at k = 0 . 

Taking into account both the quasienergy and coher¬ 
ence factors, we find that the Floquet state wavefunction 
in Eq. (II. 8 D {u^ F \e,t), v^ f \e, t)}e zE< ' F) ^ t can be well- 
captured by the approximate phase II ground state solu- 
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tion {u 2 (e), V 2 (s)}e lE2 ( e ' )t when e £_, and the excited 
state {v 2 (e), — U 2 {s)}e~ lE2( - e ' ,t when e -C £+. This is true 
even for the strong quench B. As a result, we can con¬ 
struct a full wavefunction for the phase II approximation 
as follows: 


u( II \e, t) 

/l —7 (// )(e) 

u 2 {e) 

v (e, t) 

" V 2 

M e \ 


/l + 7( /z )( e ) 

u 2 (e) 

V 2 

_-u 2 {£)_ 


7 (J/) (e) = sgn(e — £+) 7 (e), 


e -iE 2 (e)t 


(4-5) 


connected” to a special limiting case. This is the limit 
A« — 0 + for fixed A 1 ,/' 1 , which describes a quench from 
the Fermi liquid ground state perturbed by an infinites¬ 
imal seed of p + ip [H] superfluid order. In this limit 
the two isolated pairs of roots that characterize all phase 
III quenches coincide. The solution is a single soliton 
mum in the order parameter. For the p + ip case, 
this takes the form 

A(t) = y/R(t) expH0(i)], 

R(t] = _V_ (4 ' 6) 

|u| cosh [2ui (t — t 0 )] + u t ’ 


where 7 (e) is the true distribution function in phase III. 
The sgn function is necessary to “unwind” the distribu¬ 
tion function, since a phase II quench close to the II III 
border has y(e — 0 ) = — 1 0 . 


C. Phase III dynamics, avoided crossing (“Floquet 
bandgap”), and BCS instability of the normal state 

As discussed above, all phase III quenches feature a 
population imbalance in the basis of Floquet states at 
low momenta: the excited state Floquet band becomes 
occupied with unit probability in the limit e = k 2 —> 0 . 
This is depicted as the winding from —1 to +1 with de¬ 
creasing e of the phase III distribution function 7 (e) in 
Fig. [9] The mechanism for this is a combination of two 
factors. First, the Floquet band structure (quasienergy 
spectrum and Floquet state coherence factors) exhibits 
exactly one “Floquet bandgap” or avoided crossing, as 
shown in Figs. [7] and [5] In particular, the Floquet co¬ 
herence factor u^ F ^(£,t) —> 1 [w( F )(e,f)] —> 0 for both 
e —► {0,oo}. By contrast, u 2 (e) [u 2 (e)] approaches one 
[zero] for large £, and zero [one] for e = 0 ; this is the 
usual behavior for coherence factors used to parame¬ 
terize positive energy quasiparticle excitations above a 
BCS ground state. In the Floquet case, the inversion 
\u^ E ^ (e, t) | —> u 2 (e) for £ -C £ + implies the winding of 
7 (e), so as to ensure the conservation of the topological 
pseudospin winding number @. 

What is the origin of the single avoided crossing in 
the Floquet bandstructure? At first glance, this result is 
surprising, as we have obtained exact results in a system 
where the (quench-induced) drive frequency D is always 
much smaller than the bandwidth (which can be taken 
as an energy cutoff A O kp/2 0 ). The unperturbed 
(e.g., phase II approximate) spectrum folds many times 
when reduced to the first quasienergy Brillouin zone. We 
note that the periodic drive A oc (t) is not in general a 
pure harmonic [Fig. [ 6 ]. We might therefore expect the 
opening of small bandgaps whenever the folded spectrum 
approaches the zone edge, in analogy with ID bandstruc- 
tures in solid state physics. 

In fact, the single crossing can be understood as a con¬ 
sequence of the integrable BCS dynamics. The key idea 
is that any quench in phase III can be “adiabatically 


where u rj i denote real and imaginary parts of the doubly- 
degenerate isolated roots. This soliton solution describes 
the nonlinear collisionless dynamics of the Cooper pairs 
following a linear instability (exponential order param¬ 
eter growth) of the seeded Fermi liquid [lH 0 . A (t) 
grows from the seed at t -C to, reaching a maximum at 
t = to; it then decays again as t — > 00 . The peak magni¬ 
tude is determined by the post-quench coupling strength; 
in terms of one has 

u t ~ 47 rn, u\ ~ yjAnn A(/ ] , max 

valid for A[/' <C Aqcp- Here n denotes the fixed particle 
density. The decay A(t —> 00 ) = 0 occurs because the 
quenched system cannot reach the preferred ground state 
(a paired BCS superfluid) without dissipating the large 
energy injected by the quench. Ultimately, pair-breaking 
and other processes would induce thermalization, but we 
are neglecting these (see @] for a discussion of the relevant 
time scales). 

Mathematically, the soliton solution in Eq. dUD can 
be understood as arising due to the discontinuity in the 
initial spin distribution of the unperturbed Fermi step 
IT 2 I . The only effect of non-zero A™ for a general phase 
III quench is to split the pairs of isolated roots, which 
gives rise to a train of solitons in time separated by a 
finite interval T, see Eq. (TOl) . By contrast, a normal 
state initial distribution with multiple discontinuities is 
expected to seed a superposition of soliton trains with 
multiple incommensurate frequencies (l2| . 

The explicit solution of the Floquet spinors given in 
Appendix [B] combined with Eq. (14.61) implies that 

\u {F) (£,t)\ -± 1, |u (F) (£,t)| -4 0, 

in the limit that t — > 00 , wherein the soliton has de¬ 
cayed. One can also verify that | 7 (e)| = 1 for all £ when 
A ^ 5 —> 0. Thus, in the long-time limit we recover a pseu¬ 
dospin configuration consistent with a “normal state,” 
s z {e) = ±1/2. This state must possess the same par¬ 
ticle density as the initial Fermi step; moreover, given 
the connection between discontinuities and isolated roots 
discussed above, we expect only a single discontinuity in 
s z {e). We conclude that the asymptotic pseudospin tex¬ 
ture is that of the initial Fermi step. Different from a 
ground state, however, is the fact that this is encoded 
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w/£l 


FIG. 10: (color online). Time-averaged bulk rf spectra (red 
solid curves) from a topological Floquet system induced by 
the phase III quenches (a) “A” and (b) “B,” as indicated in 
the quench phase diagram Fig. [T] For comparison, the cor¬ 
responding phase II (static quasiequilibrium) approximations 
[Eq. (14.51) 1 are plotted with dashed blue curves. Unsurpris¬ 
ingly, the two curves agree well for quench A (a), which is 
close to the phase II III border. More interesting is the close 
agreement for quench B (b), which resides deep in phase III 
and is characterized by the large anharmonic oscillations in 
A«(t) shown in Fig. [ 6 ] The good agreement can be largely at¬ 
tributed to the nonequilibrium distribution functions. These 
show a population inversion in the Floquet bands at small e, 
as shown in Fig. [9] The small deviation in (b) consists of a 
sequence of evenly spaced peaks (Floquet copies 0 ) , with 
the spacing equal to the oscillation frequency f 2 . 


in the sign of the distribution function 7(e), rather than 
the coherence factors (since \u^ F \e, t)| = 1). This is 
again because the quench-induced state is at all times 
very far from the preferred ground state. The relation 
\u( f >(e, t)| = 1 for 1 t — to\ —> ±00 implies that the quasi¬ 
particle creation operator aj { ~ c^, where the latter cre¬ 
ates an elementary fermion. In other words, all fermions 
are excitations for a quench that turns on pairing interac¬ 
tions, due to the BCS instability for any kp > 0. These 
characteristics of the non-equilibrium coherence factors 
and distribution function apply to all phase III quenches. 
As can be seen for the deep phase III quench “B” shown 
in Fig. [lib), 7(e) approaches the Fermi step for small 

A$*>. 



a)/Q 



aj/Q 

FIG. 11: Time-averaged bulk rf spectra of the lower Floquet 
band associated to quenches (a) A and (b) B. These spectra 
are calculated using the same Floquet states and quasiener¬ 
gies utilized in Fig. [10] but assuming 7 (e) = —1, that is only 
the lower Floquet band is occupied. Although these only dif¬ 
fer in the distribution function, the spectra are completely 
unlike the corresponding ones in Fig. 1101 In particular, the 
case (b) corresponding to the Floquet bandstructure induced 
by quench B gives an rf spectrum from the lower Floquet band 
that exhibits no gap, in contrast to quench-induced spectrum 
in Fig. llOl bl. which takes into account the physical distribu¬ 
tion function. 


D. Phase III: RF spectroscopy 

In this subsection, we discuss the rf spectrum of the 
quench-induced Floquet system and compare it with the 
phase II approximation [Eq. (14.51) ]. Using the exact re¬ 
sults for the Floquet coherence factors and the distribu¬ 
tion function, we evaluate Eq. (EH). The time-averaged 
rf currents (red solid curves) for quenches A and B are 
plotted in Fig. [TO] together with the corresponding phase 
II approximations (blue dashed curves). As in Sec. IIII1 
here we assume that the non-pairing state is initially un¬ 
occupied for all momenta in the rf calculation, = 0. 

We expect to see two series of peaks evenly spaced 
by the quench-induced frequency of oscillation U. These 
come from processes wherein an rf photon breaks a 
ground or excited state Cooper pair, and absorbs or emits 
several oscillation quanta |18|| . No satellite peaks are vis¬ 
ible for quench A, which is quite close to the phase II 
border and is characterized by a small oscillation am¬ 
plitude [see Fig. [B], In this case the signal is domi¬ 
nated by the zeroth order term. The stronger quench 
B is deep in phase III; in this case, A oc ,(t) exhibits an¬ 
harmonic time-dependence and a large oscillation ampli¬ 
tude. Even for this quench, however, only the lowest 
order satellite peaks are visible. This is because the os- 
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FIG. 12: Real part of the bulk first order rf harmonics 

Re/ r /(l) for quenches (a) A, (b) B. 



V/Q. 



cillations in u T) (e, t ) are confined to a particular region 
near £ + < e < £_, as shown in Fig. [ 8 ] 

In Fig. UOf a), the phase II approximation (blue dashed 
curve) agrees well with the actual rf spectrum for quench 
A (red solid curve). The deviation is also relatively small 
for quench B, depicted in Fig. HHIb). We note in partic¬ 
ular that the bulk rf spectrum shows a robust gap, even 
for this strong phase III quench. We attribute this to the 
crucial role played by the distribution function. As dis¬ 
cussed above, the topology of the initial pre-quench state 
constrains q(e = 0) = 1, so that the upper Floquet band 
is occupied at low energies. The inversion of y(e) from —1 
to +1 with decreasing e effectively interchanges u <F " 1 (e, t) 
and v( F \e,t) in Eq. m- Since u^ F \e,t) behaves like 
u 2 (e) [u 2 (e)] for £ £± [£ < £±], while v' F \e,t) shows 

the opposite behavior, the combination is well-captured 
by Eq. (031) . 

To compare, we consider a Floquet system where the 
same oscillating order parameter is imposed externally. 
Unlike the quench-induced case, we assume that this sys¬ 
tem is prepared in a way so that only the lower Floquet 
band is occupied, i.e. 7 (e) = —1. The new time-averaged 
rf spectrum is shown in Fig. 111! and is found to be dra¬ 
matically different from that of quench induced asymp¬ 
totic state. In particular, for the Floquet states and 
quasienergy spectrum associated to the strong quench B, 
populating the lower Floquet band gives a bulk rf signal 
that does not exhibit a gap, Fig. [TTT bh 

In addition to the time-averaged value, the rf signal 
in the Floquet phase exhibits harmonics at the drive fre¬ 
quency. In Fig. [12] we plot the real part of the first order 
bulk rf harmonics signal, Re/ r /(l). This is defined via 


FIG. 13: (color online). Time-averaged bulk tunneling con¬ 
ductance for a topological Floquet system following quenches 
(a) A, (b) B. The red solid curve and blue dashed curve in 
this figures respectively show the phase III conductance G(V) 
and its phase II approximation. These two curves agree well 
in (a) but not in (b). The phase II approximation fails for the 
strong quench B as the distribution function does not influ¬ 
ence the tunneling signal. C.f. the rf spectra in Figs. 1101 and 

m 

Ea. (12.121) and coumuted in Ea. ED in Appendix El In 
contrast to the average spectrum, a clear sequences of 
peaks is observed for both quenches. 

E. Phase III: Tunneling 

Aside from the rf signal, we also compute the time- 
averaged tunneling conductance, assuming that phase III 
can be realized in a solid. For the quench A located in the 
vicinity of the phase II III boundary, the phase II approx¬ 
imation works well for most V, as shown in Fig. 1131 a’). 
By contrast, it fails for quench B. In particular, the gap in 
the static quasiequilibrium phase II approximation does 
not exist in the real spectrum, see Fig. fl3l b'). These re¬ 
sults should be compared to the bulk rf spectra, Fig. [TO] 
and Fig.QD the former uses the physical quench-induced 
7 (e), while the latter assumes 7 (e) = — 1. We conclude 
that the absence of a gap in the tunneling signal for 
quench B can be attributed to the fact that the latter 
is independent of 7 (e), Eq. (12.111) . 

Fig. 1141 illustrates the harmonics of the tunneling con¬ 
ductance Re G(l) [Eq. (12.121) ] for quenches A and B; these 
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FIG. 14: Real part of the bulk first order tunneling harmon¬ 
ics ReG(l) for quenches (a) A, (b) B. 

are computed via Eq. (IA2I) . While clear peaks are ob¬ 
served for the weaker quench A similar to the rf shown 
in Fig. [12l the result is much less clear for the strong 
quench B. 


F. Phase III: ARPES 


two series of Floquet copies in the spectra, and the in¬ 
tensity decreases fast as the Fourier order increases. For 
comparison, the phase II approximations associated with 
quenches A and B are respectively shown in Fig [15] 
(c) and (d). We notice the phase II approximation re¬ 
sembles the leading order copy in the ARPES signal for 
most momentum fc, but not for those around k+ = yf£+ 
[Eq. (IB20I) ] where the wavefunction exhibits the biggest 
difference from its phase II approximation [Eq. (14.51) ] , es¬ 
pecially for the relatively strong quench B. 

V. LOCAL BULK AND EDGE SPECTROSCOPY 


In this section, we look for experimental signatures of 
Majorana fermion edge modes [32|[34J in a 2D p+ip topo¬ 
logical superfluid. Specifically, we investigate the contri¬ 
bution from the edge modes to the local rf and tunneling 
signals. Our results apply equally to the BCS ground 
state, and to a quench-induced topologically non-trivial 
superfluid in phase II. The latter is characterized by a 
constant and > 0 [Eq. (11.41) ]: quenches produc¬ 
ing states of this type are indicated by the shaded portion 
of region II in Fig. |T| 

The system is assumed to occupy an infinite strip 
of width C in the x-direction, extended in the y- 
direction. Our results also apply to the half-plane 
{(x,y)\x > 0, y £ M}. At a physical edge x = XEdge, we 
impose hard wall boundary conditions. 


The time-averaged ARPES signals from the Floquet 
system induced by phase III quenches A and B are il¬ 
lustrated in Figs. [15] (a) and (b), respectively. We find 


A. BdG Hamiltonian 

The effective Bogoliubov-de Gennes mean field Hamil¬ 
tonian in this geometry is given by 



FIG. 15: (color online). The ARPES spectra for the topo¬ 
logical Floquet system induced by Phase III quenches (a) A, 
(b) B. Their phase II approximations are depicted in (c), (d) 
respectively. 


HsdG 


AS X ) c -ky(x) Kky) 


Cky ( X ) 


1 dr 


h(ky)= 2 ^ + 2 lky) . (5.1) 

A* (—i— + ik 1 iAL — tjL _J_ /; 

l dx +tK y> 2 dx 1 2 + A 


After the Bogoliubov transformation, 


' a ky ip) ' 

f 

' Uky (P, X) V ky (p, X) ' 

' C k y(X ) ' 

a ~ky(P). 

~Jx 

V-ky(P,X) U*_ k J P ,x)_ 

_ C -ky( X )_ 


where the amplitudes satisfy 


h*(ky) 


Uk v {p,x) 
Vky ( P , X) 


= E ky ip) 


Uky{P,X) 

Vky(P,X) 


with Ek y ( p ) > 0, the Hamiltonian reduces to 


Hb6G =2 


5 Z Ek v(P ) 


^y p 


a ky(P ) a ky(P ) 

- a-ky{p ) al k (p) 


(5.2) 


(5.3) 


(5.4) 
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Here a ky ( p ) and a\ ( p ) are the annihilation and creation 
operators associated with the quasiparticle excitation la¬ 
beled by p and transverse momentum k y . The corre¬ 
sponding quasiparticle energy is denoted as E ky ( p ). 
runs over all the eigenstates with nonnegative eigenval¬ 
ues, including edge states (if present). 

Atoms in the non-pairing hyperfine state |2) in the rf 
experiment are described by the Hamiltonian 


H { 0 d) = 


4„ {x)hd{k y )dk y Or), 


1 d 2 fc 2 

hd{ky)= ~2dx2 + il 


(5.5) 


The local current I(jo,t) therefore is composed of two 
components from the commutator with Ht and , 
respectively: 


I( ro) =It{ r 0 ) +Id{ro), 
I T ( r 0 ) =i[H T ,n d ( r 0 )] , 


Id{ r 0 ) =i 


H (d \n d (r 0 ) 


(5.9) 


In the interaction picture, the current expectation value 
is given by 


which can be reduced to 

H o d) = \ E [^ y (k)+S2,i] b{ y (k)b ky (k), ( 5 . 6) 

Jk v k 

where b ky (k) annihilates an atom in state | 2 ), labeled by 
k and with energy £ ky (k) + £ 2 ,i- 

h y (k) = 

h* d {ky) w{k , x) = [£ fey (k) + £ 2 , 1 ] w{k , x). 

Here£fc y (fc) = {k 2 +k y )/2—p,^. The wavefunction w(k, x) 
for the non-pairing state is a real sinusoid that vanishes 
at the edges of the cloud. In what follows, for clarity we 
use p and k to indicate the excitation in the paired and 
unpaired atomic levels, respectively. 


/ 


w(k,x) d ky {x), 


(5.7) 


B. Local rf current 

From the Heisenberg EOM, the local current at ro 
equals 

(5g) 

H = H BdG + H { 0 d) + H t . 


{U\t) I{t) U{t)) 0 = h+I 2 + I 3 + O (T 3 ) , (5.10) 

where 


h(t)=-if dt'{[lT{i),H T {t')}) o, (5.11a) 

J —OO 

hit) = f dti f dt 2 (HT(ti)I d (t)HT{t 2 ))o, (5.11b) 

J— OO J— OO 

{HT(tl)HT{t 2 )I d (t))o 
+ (Id{t)HT{tl)HT(t2))o_ 
(5.11c) 



In these equations, (...)o indicates expectation with re¬ 
spect to the state without the rf perturbation. In the 
homogeneous case, only Ii(t) takes a nonzero value. 

We will assume that the initial density of non-pairing 
state 12 ) atoms is equal to zero; it is straightforward to 
treat the more general case, but the result is rather cum¬ 
bersome. In this case the total current is positive definite, 
while h{t) in Eq. (15. 1 lcl) vanishes exactly. For the finite 
sample geometry of interest here, both h{t) and I 2 (t) 
are comparable and must be evaluated, e.g. Ii(t) can 
take negative values for a system with an edge. We find 
the following result for the time-averaged local current at 
x = Xq: 


J 


I(x 0 ,u) = 2ttT 2 / y 


2 t 1 +7 k y {p)\ \w{k,x 0 )| 2 

J u ky (p,x) w*(k,x) 

6 [u + E ky (p) -£ ky (k)] | 

| +2 [!-7 k y ip)\ \w{k,x 0 )| 2 

J v ky (p,x) w*(k,x ) 

2 

5 [w - E ky (p) - ik)] 

j 


(5.12) 


In this equation, 7 k y [p) = 2 (al y {p) a ky {p)) - 1 encodes 
the occupation of the states [induced by the quench; c.f. 
Eq. e^l- 

In the derivation of the local current, we did not use the 


explicit form of Uk y {p,x ), Vk y (p,x), or w{k,x), but only 
the assumption that each is either purely real or imag¬ 
inary. This assumption should be true for any realistic 
boundary conditions given the form of paired and non- 
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pairing Hamiltonians, Eqs (E3D and (1531) . Eq. (15T21) 
therefore applies to any semi-infinite or infinite system 
governed by Eq. 

Eq. ( 15 . 121 ) takes the form of Fermi’s Golden rule. The 
first (second) term is due to the process where the excited 
(ground) state Cooper pair with energy E ky ( p) [—E ky (p)] 
absorbs an rf photon, producing a state | 2 ) atom with en¬ 
ergy £, ky (k) and an unpaired state | 1 ) atom with energy 
zero. The factor \ [l + 7 k y (p)\ {\ [l - 7 k v (j>)]) encodes 
the probability that the quasiparticle state {p,k y } is ini¬ 
tially occupied (unoccupied) before the application of rf 
radiation. The factor | w(k, £o)| in each term is the prob¬ 
ability that an atom exists at x = xq in the non-pairing 
| 2 ) state with momenta {fc, k y }. 

Although Eq. (I5.12[) is very similar to the homogeneous 
result [e.g. Eq. (12.91) with only Mo.kjfio.k nonzero], it has 
very different implications for bulk versus edge states. 
The transition rates are proportional to the squared- 
overlaps 


T n-ky (p, k) = 


Tv,k y (p,k) = 


J u ky (p,x)w*(k,x) 
/ v ky (p,x)w*{k,x) 

J X 


(5.13) 


When coherence factors u ky (p, x) and v ky (p, x) describe 
extended bulk states, then the product of these with 
extended (box) wavefunction w*(k,x) is of order 1/C, 
where C is the extension of the system in the x-direction. 
For p = k, the integral over x gives a factor of C, so that 
T u ,v-k y are of order unity. When u ky (p, x) and v ky (p, x) 
describe a spatially localized edge mode, these are a fac¬ 
tor of 1/C smaller. We conclude that the contribution 
of the edge states is suppressed by a factor of the lin¬ 
ear system size. Moreover, because k y runs over a finite 
range for the edge modes, the edge signal is delocalized 
throughout the bulk of the sample (since it obtains from 
integrating \w(k, xq )\ 2 over a finite range of k for which 
energy conservation is satisfied). 

Thus, local rf spectroscopy of this type is not well- 
suited to detect Majorana edge modes. 


C. Tunneling 


by 


By contrast, the local tunneling conductance is given 


G(x 0 ,V) 


= 2ttT 2 v 0 / 


\uk y {p,x 0 )\~ 6 V + E ky (p) 


p +\Vk y (p,x 0 )| S V — E ky (p) 


(5.14) 


This is independent of the distribution function (qucnch- 
induced occupancy), and gives equal weight to discrete 
bound states (edge modes) and continuum bulk modes. 
The contribution of each of the latter is suppressed by a 
factor of 1/C, but this is compensated by the summation 
over p. 

When the metal tip is placed deep in the bulk, the tun¬ 
neling spectrum resembles that of the homogeneous case. 
Near the boundary, signal within the gap is contributed 
only by the edge modes. 


VI. CONCLUSION 


In this paper, we have obtained the rf and tunneling 
spectra of a quench induced out-of-equilibrium steady 
state. In particular, we focused on the Floquet system 
with time periodic order parameter induced by a phase 
III quench. One important aspect of our result is the 
connection between these experimental observables and 
the distribution function. We found that the distribution 
function plays an essential role in the bulk rf spectrum. 
We demonstrated that the reason the time averaged rf 
(but not tunneling) spectra are in good agreement with 
static quasiequilibrium (phase II) approximations lies in 
the distribution function, which is forced to exhibit a 
population inversion in the basis of Floquet states due to 
the conservation of the topological pseudospin winding 
number. This crucial information is missing in the tun¬ 
neling spectrum, and leads to the disappearance of the 
gap in tunneling spectra for strong phase III quenches 
characterized by large amplitude oscillations in the or¬ 
der parameter. Finally, we showed that local rf is not a 
good method to detect Majorana edge states, due to the 
non-local character of the radiation-induced transitions 
to non-pairing states. 
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Appendix A: RF and tunneling amplitude 
harmonics 


In the superfluid Floquet phase, the rf current given by 
Eq. (12.71) exhibits modulations at harmonics of the drive 
frequency Q. These are encoded in 
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Mp) = ^£ 


1T 2 


E 


(1 - 7k)(l - ft n ) ~ (! + 7k)/k° (v^ M v n+pM + v*_ pM v nM )5(uj - Ck - £'k F) + nfi ) 
+ (1 + 7k)(l - /k d) ) - (1 - 7k)/k d) (^ )k «n+p,k + w*_p, k u„,k)^(w - Ck + -Bk F) - 
(1 - 7k)(1 - /k d) ) - (1 + 7k)/k rf) ] («n,kUn+p,k - «n-p,kVn,k)- 


W — £k — -E^ + 


(1 + 7k)(l - /k d) ) - (1 - 7k)/k^ K,k“™+P,k - <-p,k«n,k) 


(d) 


w - £k + -Ek F) - - 


y (ai) 


An analogous expression for harmonics of the tunneling conductance is given by 

G(p) =7rT 2 r'o E [(^n.k^n+P.k + v*-p,kVn.k)S{V - + nfl) + (u* k“n+p,k + u^_ pM u n ,k)S{V + eI F) - nil) 

1 , , 1 


- iT 2 v o ^2 

n,k 


{^n,k^n+p.k V n -p.kVn,k) 


V - eI F) + nn 


(u n k'Un+p,k U n —p. k^n.k) 


V + eI F) - nn 


(A2) 


We find that u n> k 7 Vn.k are real, so the second term in 
Eqs. (EH) and (IA2I) can be ignored if only the real parts 
of I r f(p ) and G(p) are needed. 

Appendix B: Floquet states and occupations via 
integrability: Explicit solution in Phase III 

Through self-consistent mean field theory and the use 
of the Lax construction, the order parameter A(t) as a 
function of time is determined (A = The am¬ 

plitude (VR) and the argument (_—<P) of the complex or¬ 
der parameter follow the EOM [9j 

R 2 = (R+ - R)(R - R_)(R + R+)(R + R_), (Bl) 


^=^R + 2m-|. (B2) 

where 

R± = (Ar ± A 2 ) 2 , 

/ /- /- \ 2 (B3) 

R± = f Y 2/^1 ^ A 2 ± y 2fi 2 - A| J , 

m = Emi + M2 - A? - A|) , 

2 (B4) 

= - (A 2 - A 2 ) (A 2 - A 2 - 2yur + 2 M2 ). 

The parameters Ai 2 and are determined via 

Eq. (14.11) by the two pairs of isolated roots that char¬ 
acterize a phase III quench. For a given quench specified 
by the coordinates {A™, A ( / ] }, these roots solve a cer¬ 
tain transcendental equation [9f. Results can be obtained 


numerically for a given fixed particle density n and ul¬ 
traviolet energy cutoff A. 

The solution to Eq. (IB2I) is a combination of a 
periodic part 4>(t) and linear part 2 n^t. Here $(t) shares 
the same period T with R (f), and is defined as 

hi dt *=2 hi 

(B5) 


a. Lax reduced solution 

In the following, a “Lax reduced” solution Q to the 
BCS spin dynamics is utilized to find the explicit solution 
to our problem, i.e. the asymptotic steady state following 
a sudden quench of coupling strength. 

We consider a reduced solution with the same order pa¬ 
rameter (isolated roots) as a particular phase III quench, 
but with different initial conditions. The existence of 
such solutions, along with an explicit prescription for con¬ 
structing them, was detailed in Secs. Ill C and IV B of 
0, based upon previous s-wave work [3^J . In terms of 
Anderson pseudospins, the reduced solution can be writ¬ 
ten as a product of a sign function £(e) and a pseudospin 
function so(e), 

sr(s) =C(e)so(e), (B6) 

where £(e) = ±1. The single particle levels are labeled 
via the momentum-squared e = k 2 . To qualify as a so¬ 
lution, £(e) will exhibit a discontinuous jump at some e 
from +1 to —1; this is closely related to the soliton so¬ 
lution discussed above in Sec. IIV Cl The function s 0 (e) 
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depends only on the order parameter, and is given by 
sg(e) =a(e)R + 6(e), 


s o ( £ ) = 


Q(g) 

2y/eR L 


—zR + R 2 + (4m — 2e)R + 2 ip 


where a(e) and 6(e) are functions of e, 

«( £ ) = 


b(e) = - 


4E 1 (e)E 2 (ey 


(B7) 


(B8) 


m{e)E 2 {£) 


Here E a (e) = i/(e/2 - fj, a ) 2 + A 2 e, a = 1, 2 is a p + ip 
quasiparticle energy [Eq. (11.31) ]. 

The order parameter is defined self consistently by 


A(f) = -Gj2VeiSi ■ 


(B9) 


Substituting in this equation from the reduced solu¬ 
tion in Eq. (IB6I) and using Eq. (IB71) . we find that the 
L.H.S. reduces to 


- G(£ C^) ~' R+R2+ ^ mR + 2V, e~^ 
i VR 

+ G(^Ci£iai)VRe"**. 

i 

This gives \/Re~ 1 ^ if and only if 

'y \ Ci a i = 0) 'y \ c 


(BIO) 


iSidi - G , 


(BH) 


which serve as the constraints for the Lax reduced solu¬ 
tion. Under such constraints, we find the total z-spin J 
(related to the fixed number of particles N) and the total 
energy E are indeed conserved quantities since 


j=T,Ci4,i= Y.tibi, 

i i 

_^ p _^ 

E = ^ Ci £ i s 0,i — ~q = 


(B12) 


Another important quantity is the Lax norm Q, which 
generates all the integrals of motion. We find that the 
Lax norm of the reduced solution is given by 


L 2(w) = ( 


Qi~i 


2y / Q 4 (e i ) £i - u 


Q 4 (u). (B13) 


The {Ci} can in principle be determined by equating 
Eq. (IB13I) to its initial value. 

To find the spinor of So parameterized as 


(B14) 


u 0 {e,t ) 


'|uo(M)|e~ w “ {e,t) ' 

v 0 (e,t) 


Juo(e,t)|e- !;e ”( £ ’ t )_ 


we use 


s 0 (e,t) =Uo{e,t)v 0 (e,t) 

s o( £ ’ l )=\ (lw( £ , t) | 2 - | u 0 {e, f)| 2 ). 

In the rotating frame s~ —> s~e 1 ^, this yields 


(B15) 


M £ , t)\=\j~- [a(e)R(f) + 6(e)], 


(B16) 


l^o( £ , *)| = - + [a(e)R{t) + b(e)], 


8 u (e, t) — 8 v (e, t) = arg zR + R 2 + 4mR + 2ip — 2eR^ . 

(B17) 

Inserting Eq. (IB16I) into Eq. (11.71) . we get 

n £ , , a(R 2 +4mR + 2^-2eR) 

9u -~2 + 2M + 1 — 2 (aR + 6) ’ (B18) 


• el# a (R 2 + 4mR + 2-0 — 2eR) 
v 2 2 dt + 1 + 2 (aR + 6) ' ^ ^ 


We define and e- by 


Rj_ + 4mR-{- + 2-0 

£± ~ 2r; ■ 


(B20) 


Eq. (IB18I) [Eq. (IB19I) ] is valid as long as e ^ e+ or t ^ 0 
(e y e_ or t ^ -j). 0„(e + ,O) and 0„(e_, -j) are unimpor¬ 
tant since |zz(e+,0)| = 0 and |v(e_, ^)| = 0. 

We next switch from the frame rotating with phase 
(f>(t) = 2/z^t + <I>(f) to that rotating only with its linear 
piece 2/z^t. This is equivalent to replacing d<j>/dt —> 2/z^ 
in Eqs. (IB18I) and (IB 191) . The resulting equations depend 
on t only through R, which implies that 9 U (6 V ) should 
be a sum of a time-periodic function and a linear piece 
— E u t (—E v t). The intercept E u (E v ) is determined by 
integrating Eq. (IB18I) [Eq. (IB19I) ] over one period, 


1 


E u {e) = - — dt6 u {e,t) 


T 


dt /_(e, t), 


(B21) 


E v (e) = - — dt 9 v (e, t) 

1 Jo 

= dtf+{e,t). 


Here 


cl ^ _ a( £ ) [R 2 (i) + 4mR(t) + 2?/)-2eR(t)] 
■ f±{£ ’ t >~ 1 ± 2 [a(e)R(i) + 6(e)] ' 1 ’ 
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We find that E v (e) — E u (e ) = —no(e)fi from 
Eq. (IB17I) . where no(e) = 9(e — e + )9(e_ — e). 
As a result, {u 0 (e,t), vo(s,t)} is a Flo- 
quet state and can be parameterized as 
{|uo(e, t)\e~ ie>u< ' e,t \ \vo{s,t)\e~' iev( ' e ^}e +lEo ^ t . Here, 
the quasi-energy Eq = E u (mod Q), Q u (t + T) = Q u (t) 
(mod 2n), and 0„(f + T) = Q v (t) (mod 27r). 0„,0„ 
can be obtained by integrating Eqs. (IB18II and (IB19I) . 
Assuming R(t = 0) = R + and $(f = 0) = 0, we deduce 
the integration constant from Eq. (EUD, 


0 „(e, 0 ) - 0 „(e, 0 ) = 7 r 0 (e - e+), 
T T 

-) - j) =7 T0 (£-£-), 


(B23) 


Then we obtain 


&u{£>t) — J dt'f-(e,t')+ — — + /Too + E 0 (e) 

Q v (e,t) = - tt9 (e - e + ) + f dt'f + (e,t') 

Jo 


12 


— Moo + Eo(e) 


(B24) 


In the scheme where we take the quasienergy Eq = E u , 
we have Q u (e,t + T) = Q u (s,t), 0„(e, t + T) = ©„(e,t) + 
27rno(e), and Eq. (IB24I) simplifies to 


0«(e,i) = 

@v(e,t) = 


[ dt '~y» dt1 


f-(£,?), 


dt,' — 


T 


dt' 


+ 2nd ( e — e + ) 6 (e_ — e) 


f + (E,t')-n0 (s-e + ) 
t 


T 


(B25) 


Combining Eqs. (IB 161) and (IB24I) . we obtain the complete 

expression for [uo(£,t) vo(e,<)] T . Then the spinor of the 
Lax reduced solution can be written as 

K(M)|e- i0 “ (e ’ tr 


u(e, t) 

/ 1 +C(e) 

v(£,t) 

V 2 


o +iE 0 (e)t 


l-CC) 

2 


' \v 0 {E,t)\e ie ^ 
-|u 0 (e,t)|e ie “^. 


0 —iEo(e)t+iV(e) 


(B26) 


Here T(e) is some time-independent phase. 


The Lax reduced solution in Eq. (IB26I) is an extreme ex¬ 
ample where 7 = — C only takes values equal to ± 1 . In 
order to solve the self-consistent conditions in Eq. (IB11I) , 
((e) must exhibit a discontinuity (“Fermi step”). The re¬ 
duced solution therefore cannot apply to a quench from 
a BCS initial state, since the latter initially has a smooth 
pseudospin texture in momentum space, and the BCS dy¬ 
namics do not change this. The solution for the quench 
involves replacing —((e) with a smooth 7 (e) that never¬ 
theless winds from +1 to —1 with increasing e. 

The conservation of energy E, ^-component angular 
momentum J (particle number), and the Lax norm al¬ 
lows us to determine the 7 j for the post-quench state. 
Substituting Eq. (IB27I) into Eq. (IB 151) . we have 


s i (*) = - HU* 0 A v 0,i 




( ■_^ iu *2. e -2iE i t + e -iT iv 2 ie 2iE 4 t 




(-« 


1 

—‘liEit 


e u 0 ,i V 0,i e 


(B28) 


These can be considered as a combination of so,i(t) 
and an oscillating spin with energ y-depen dent frequency 
s 0 sc,i(i)) weighted by —% and ^ 1 — 7 f separately, i.e. 
Si(t) = -7 iSo,i + \/l - 7 i Sosc ,i- As t -*■ 00 , through re¬ 
peated integration-by-parts or the saddle point approx¬ 
imation, we find that the contribution of s osc ,i to E, J , 
and the Lax norm vanishes. Therefore, all of the con¬ 
served quantities can be obtained by replacing ( with 
—7 in Eqs. (IB11I) (IB13I) . which leads to 


XI7 i<h = 0, 

i 

X ^ = - J > 

i 


X Ji £ i a i 

i 



x^=- ^ 


and 


L2(u) 


l 7 j £ j _ 1 

2^Qj£~)£j - u 


2 

Qa(u) . 


(B29) 


(B30) 


b. Initial condition 


We assume that the general solution to the self- 
consistent BdG equations following a quench takes a sim¬ 
ilar form, 


u(e, t) 

1 

T 

y(e,t) 

V 2 


. 


\u 0 (s,t)\e- i@ 

|u o (e,t)|e i0 ^ E ’ t ) 

-|uo(M)|e i0u(eit) . 


e +iE 0 (e)t 


e -iB 0 (e)t-|-ir(e) 

(B27) 


We now let u in Eq. (IB30I) approach the positive real axis 
from above and below, u —> u ± ir], u G R + and rj — > 0 . 
From the residue theorem, 


L 2 (u ± ir]) = v 2 


P r de ■ 7(E)e 


1 


/0 2yJ Q 4 (e) £ ~ u 

|2 (B31) 


±17T 


"f(u)l 


2y/Ol(u) 


£U(w), 


where ^ = £ 2 / 87 t is the bare density of states, and £ 
denotes the linear system size. This is equal to the Lax 
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norm of the pre-quench state, 


order parameter is given by 


L 2 (u ± irj) 



which gives 


\/R (t) = A 2 (const.), 
<f>(t) = 2 n 2 t. 


(B34) 


l7(e)l 


\y/T^e)-VUe)\ 

27re 


The Lax reduced solution on the border can be obtained 
' ^ ' from Eqs. (1B34I) . (1B7I) and (1B8I) . 


This shares the same form as the pseudospin distribution 
function in phases I and II. Here, /+(e) is a certain func¬ 
tion that depends upon the parameters of the pre-quench 
state, as well as the strength and direction of the quench 
(i.e., the quench coordinates {A™, A^/ 5 }), while /-(e) is 
its complex conjugate; an explicit expression appears in 
[§|. The sign of 7 is determined by enforcing continuity 
whenever 7 (e) — > 0 with a non-zero slope [9|. We find 
that 7 (e) winds exactly once from +1 at small e to — 1 
at large e for all quenches in phase III; see e.g. Fig. [51 
For all phase III quenches that we have investigated in¬ 
cluding “A” and “B” indicated in Fig. [T] we have used 
the result for 7 (e) and verified that all four relations in 
Eq. (IB291) hold. 

We find that F(e) in Eq. (IB27I) must be 0 or tt due 
to the invariance of the pseudospin equations of motion 
under the effective time-reversal transformation [ 9 } 

s (e, t) —^ s (e, —t), 

s ± (e,t) —t s =F (e,-t), (B33) 

A (t) -7- A*(—t). 


s z (e) 

s“(e) 


« E )sgn( £ -2 w )lLg> 

C(e) sgn(e - 2 m) ^ e~ 1 2 3 4 5 6 7 8 9 10 ^ 2 * 
2 E 2 {£) 


(B35) 


These are the pseudospin configurations for ground (ex¬ 
cited) state Cooper pairs when ((e) sgn(e—2/ri) = 1(— 1). 

In addition, following all the equations for the spinor 
in phase III, we find 


E 0 (e) = sgn(e - 2m)E 2 (e), 


1 / I /l sgn(e- 2 /xi)(f - m) 

K(e)l= V2+— im — 1 


M^)l = 


'1 sgn (e — 2 /ri)(| — fi 2 ) 


(B36) 


, 2 2 E 2 (e) 

0 u (e) = 0, ® v (e) = -tt9(s - 2m), 


c. Phase II, III border 

On the border of phases II and III, Ai approaches 0, 
and the time evolution of the modulus and phase of the 


which is consistent with pseudospin configuration 
[Eq. (IB35I) ]. Since e + = e_ = 2/x 1; E u = E v holds for 
any value of e. It is obvious that Eq. (IB27I) will take 
the form of the phase II wavefunction after applying the 
transformation 7 — > sgn(e — 2m)l- 
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